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A MATHEMATICAL MODEL OF THE CH-53 HELICOPTER 


William R. Sturgeon and Janes D. Phillips 
Ames Research Center 


SUMMARY 


A mathematical model suitable for veal-time simulation of the CH-53 heli- 
copter is presented. This model, which is based on modified nonlinear classi- 
cal rotor theory and nonlinear fuselage aerodynamics, will be used to support 
terminal -area guidance and navigation studies on a fixed-base simulator. 
Validation is achieved by compering the model response with that of a similar 
aircraft and by a qualitative comparison of the handling characteristics made 
by experienced pilots. 


INTRODUCTION 


Terminal-area guidance and navigation helicopter research is to be con- 
ducted at Ames Research Center. Prior to actual flight tests, advanced con- 
cepts and procedures will be evaluated using a piloted flight simulator. This 
simulator facility consists of a "fixed-base” cockpit, configured to that of 
the CH-53 (fig. 1), and a Sigma 9 digital computer. Operation of this simu- 
lator requires the use of a CH-53 mathematical model that can operate in real 
time on the Sigma 9 host computer. 

Helicopter models range in complexity from linear models, which are valid 
near one particular flight condition, to nonlinear blade-element models which 
account for complex rotor flow conditions and are used over the entire flight 
regime. A model of Intermediate complexity, which meets simulation require- 
ments for terminal-area guidance and navigation studies, is based on quasi- 
static rotor representations. A CH-53 model of this latter type is presented. 

The help of the following persons in obtaining this mathematical model is 
acknowledged: Dean E. Ccoper, Thomas H. Lawrence, and Phil Gold of Sikorsky 

Aircraft Division of United Technologies, Stratford, Connecticut; and J. D. 
Shaughnessy of Langley Research Center. The model was programmed on the 
Sigma 9 computer by Boris Voh of Computer Science Corporation. Validation was 
performed with the help of George Tucker and Ron Gerdes oi Ames Research Center. 


MATHEMATICAL MODEL 


The helicopter mathematical model is defined in terms of submodels of the 
fuselage aerodynamics, rotor systems, engine and governor, and control system. 
The relative relationship of these submodels is discussed in the section 
entitled "General Model Description" which precedes detailed descriptions of 



each submodel. The submodels are defined In terms of forces, moments, and 
motion expressed In the following coordinate systems which are used In the 
development of the mathematical model (fig. 2 ). 


Coordinate Systems 

1. Earth axes, subscript e: Origin fixed on the Earth's surface, x^ 

axis pointing north, y^ pointing east (fig. 2 (a)). 

2. Helicopter body axes, subscript h: Origin at the center of gravity 
(c.g.). X{^ axis forward in the plane of symmetry and parallel to the water- 
line. axis down in the plane of symmetry (fig. 2 (a)). 

3. Shaft axes, subscript s: Origin at the rotor hub, Xg axis rotated 

through the longitudinal shaft tilt angle 6 ^ about the yj^ axis, y^ axis 
rotated through the lateral shaft tilt angle about the Xg axis, Zg 

axis coincident with the rotor shaft (fig. 2(b)). This applies to both the 
main and tail rotors. 

4. Control axes, subscript c: Origin at the rotor hub, Zf. axis directed 

toward the fuselage along the axis of no-feathering (an axis perpendicular to 
the swashplate) , x^. axis points into the relative wind so that the yj. com- 
ponent of the relative wind is zero (fig. 2(c)). This applies to both the 
main and tail rotors. 

5. Wind-tunnel axes, subscript wt: Origin at the wind-tunnel mounting 

point, x^,^ axis pointing into the relative wind, z^^^ down and perpendicular 
to the relative wind. 


General Model Description 

The helicopter model is defined in terms of th« following submodels: 

1. Fuselage aerodynamics model: The fuselage aerodynamics model defines 

nonlinear lift, drag, and side forces as well as pitching, rolling, and yawing 
moments in terms of a wide range of fuselage angles of attack and sideslip, 
rotor downwash, body angular velocity, and dynamic pressure. 

2. Rotor model: nonlinear models for the main rotor and tail rotor 

define thrust, drag, and side forces as well as hub force and moments repre- 
sentative of articulated rotors over a wide range of airspeeds through hover 
to rearward and sideward flight. The rotor models account for variable inflow 
velocity, variable rotor speed, blade twist, tip loss, blade coning, blade 
flapping, flapping-hinge offset, and tail-rotor 63 hinge. 

3. Engine model: An engine and governor model adapted from a heavy lift 

helicopter simulation provides a realistic time delay between aerodynamic rotor 
torque and the resulting reaction torque applied to the fuselage. The model 
includes the effects of gas turbine, power turbine, rotor inertia, and shaft 
compliance. 
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4. Control syst«D model: The helicopter control system model converts 

pilot's cyclic control, collective control, and pedal Inputs Into main and 
tall rotor cyclic and collective pitch Inputs. An automatic flight control 
system (AFCS) Is Included which provides helicopter rate and attitude stabili- 
sation xtt roll, pitch, and yaw. 

Wind and gust Inputs to the helicopter model are provided, as well as t^ 
pilot control Inputs. All forces and moments acting on the helicopter are 
outputs of the fuselage aerodynamics and the rotor systems submodels. Fuse...:ige 
forces and moments are calculated In wind-tunnel axes and transformed to body 
axes. Rotor forces are calculated In control axes and transformed to body 
axes, and the rotor moments are calculated In shaft; axes and transfcmed to 
body axes. 

The equations of motion use the total forces and moments, In body axes, 
to calculate the translational and angular body axis accelerations. The trans- 
lational acceleration is Integrated to give body inertial velocity which Is 
transformed to Earth axes and Integrated to obtain helicopter position. The 
angular acceleration is Integrated to give body angular velocity, which is 
transformed to Euler angular velocity and integrated to obtain helicopter 
attitude. 

The relative relationship of the submodels Is shown in Figure 3(a), and 
the Inputs and outputs of each submodel are shown In figures 3(b) through 3(g). 
The model parameters are given In Table 1. 


Fuselage Aerodynamics 

The fuselage aerodynamic data are given in both equation and tabular 
form. The forces and moments are given in wind-tunnel axes in terms of local 
angle of attack, local angle of incidence at the tall, sideslip angle, body 
angular velocity, and dynamic pressure. 

Airspeed In body axes- The hel’ copter airspeed Is expressed In terms of 
Its Inertial velocity and the wind velocity as 
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The free-stream angle of attack and sideslip angle are defined 


■ tan 


•‘(S) 


-IT < a^ S w 


( 1 ) 


( 2 ) 


and 




( 3 ) 


3 





respectively! where 


V„_ ^ + V' + ^ 


and Che free-strean dynamic pressure Is 
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(5) 


Main rotor downwash effect* The effect of the main rotor downwash on the 
local angle of attack is accounted for by the rotor downwash factor (ref. 1): 
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Cimi \n* rotor parameters defined in the following section. The 

fuselage local angle of attack is 

“ft ■ “f * 'mr'kt • " S Of, S * <’) 

and the local incidence at the tail is 

^t " ^to ’ ^®kt ■ ®kf^®mr 

where ej^^ and e]^£ are empirical constants. The wind-tunnel yaw angle is 

♦»c ■ -«f <’> 


Fuselage forces and moments in wind tunnel axes- The fuselage forces and 
moments in wind-tunnel axes are provided through the wind-tunnel data given in 
figures 4 through 13. These curves are entered with the fuselage local angle 
of attack (eq. (?))• local incidence at the tail (eq. (8))t wind-tunnel yaw 
angle (eq. (9)), and dynamic pressure (eq. (5)) as determined from the equa- 
tions noted. 


Since the wind-tunnel data do not cover the full range of angle of attack 
and sideslip, it is assumed that force and moment coefficients remain constant 
beyond the limits of these angles for which data are given. This assumption 
should not significantly degrade the model performance, for large values of 
the above angles generally occur at low airspeeds where fuselage forces and 
moments are relatively small. 


The fuselage forces and moments are determined as follows: 



f 

! 


4 



where 


AD 


Also, 


- 300 sin (ft^> 

- 27.9 sin (m^) 


>vt- v-r""T"j’ 

(¥)< 

/^^iwt ^ ^^2W1;^ _ 

/ '«iwt ^ f|wt j . 


*wt 

Mwt - 
\t - 


i¥h 


Transformation of fuselage aerodynamic forces to body axes- The fuselage 
aerodynamic forces are transformed from wind tunnel to body axes. 
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Transformation of fuselage aerodynamic moments to body axes- The total 
fuselage aerodynamic moments Include the basic wind-tunnel moments, additional 
moments due to the wind-tunnel mounting point being offset from the c.g. , 
damping due to angular velocity, and rotor downwash on the tail. T.n body axes 
these moments are 
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and is the main rotor thrust. Both the damping equation a..d rotor down- 
wash moment coefficient were obtained from an unpublished Sikorsky Aircraft 
report , 


The Inputs and outputs of the fuselage aerodynamic model are shown In 
figure 3(b). 


Rotor Models 

The rotor forces and moments are calculated using nonlinear classical 
rotor theory, specifically a modified Bailey representation used In refer- 
ence 1 and discussed In references 2 through 5. Important aspects of this 
rotor model are 

1. Uniform Inflow over the rotor disk Is assumed 

2. Compressibility anl b-tall effects are neglected 

3. Lagging motion of the rotor blades Is neglected 

4. Only first harmonic motion of the rotor blades Is considered 

5. The blade coning and flapping angles are assximed quasl-statlc 

This relatively simple rotor model Is used to facilitate Its use In a 
real-time simulation. This model Is valid for forward flight to about 
120 knots, hover, rearward and sideward flight to about 20 knots, auto rota- 
tions, and large-angle maneuvers. Although the model Is adequate for guidance 
and navigation studies at airspeeds greater than 120 knots, its handling 
characteristics fidelity Is degraded due to the Increasing effects of compres- 
sibility and the reverse flow region. 

The following discussion applies to both the main and tall rotors, except 
where noted. Specific application to either the main or tall rotor Is Indi- 
cated by the subscripts m or t respectively. 
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Airspeed of rotor hub In control axea~ The total airspeed of the rotor 
hubt in control axes* is required for calculation of the rotor forces and 
moments. This airspeed is initially determined in shaft axes, using the heli- 
copter airspeed and angular velocity, and then transformed to control axes 


where 
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is defined by equation (1). The airspeed at the hub is transferred into con- 
trol axes using the rotor orientation angle 
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which is obtained using the definition of control axes, that is 

Vc 5 0 

and using small angle approximations for the main rotor cyclic control Inputs 
(swashplate angles), ky and (see fig. 2(c)), 
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Hote chat the tail rotor does not have cyclic controls » and therefore the 
corresponding k{ and Bx ere zero. 


Rotor tip speed and induced flow ratios- The rotor forces and moments are 
functions of the rotor tip speed and induced inflow ratios* which are defined 
in terms of Che hub airspeed in control axes* as 

«c 

“ * mt <"> 


and 


X 



(18) 


respectively. The induced inflow ratio v is obtained by filtering the 
8teady*-state value of v. The resulting differential equation is 



•) 


(19) 


The time constant tv is included Co account for the lag associated with 
changes in rotor inflow. Note that the thrust coefficient Cf» defined below, 
and inflow ratio X are functions of v, so that equation (19) is a first- 
order nonlinear differential equation. 


Rotor thrust and coning angle in control axes- The rotor thrust, in con- 
trol axes, and the coning angle a^ arc calculated to Che third power of the 
tip-speed ratl> according to Che following relations obtained from refer- 
ences 1 and A: 


T - bcRp(RO) 




and 
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C, 
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and 60 is the effective blade pitch angle at the root (collective pitch 
angle), and 6^ is the blade twist. Note that a term involving the blade 
mass moment in equation (22) of reference A has been neglected in equation (21) 
above, for it contributes less than O.S” and is essentially constant (ref. 1). 
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Rotor flapping angles In control axea - The calculation of the rotor flap- 
ping angles requires the fuselage angular 'velocity expressed in control axes: 
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where the transformation matrices are defined by equations (14) and (16). 

The flapping angles and (fig. 2(d)) are calculated in control axes 
according to formulas obtained from reference 1. 
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For a blade with linear twist and constant chord, it can >^e shown that replac- 
ing the 8 q appearing in the references with the pitch at 3/4 radius t'^ 
and dropping will have a negligible effect on the overall solution 
(ref. 1). The pitch at 3/4 radius is 

®o.7s " ®o + 0-75 01 (26) 


Rotor drag force in control axes- The downwind component of the rotor 
force, in control axes. Is 


H • Ta' 


(27) 


where the small angle a* is a function of the useful and Induced rotor drag- 
lift power and inflow (ref. 1), but behaves similarly to the longitudinal flap- 
ping angle ai« An expression for a', which Includes the effects of fuselage 
angular velocity (ref. 1), is 
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is 


Rotor side force in control axes- The rotor side force, in control axes 
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';nli 'iquatlon, derived from equation (3) in reference 3» neglects angular 
veto ity terms and uses the previous assumption involving the pitch at 3/4 
r ad i us ^ 60 75 * 


Transformation of rotor forces to body axes- The rotor forces in control 
axes, given by equations (20)* (27), and (.29), are transformed to body axes. 
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lihere the transformation matrices are defined by equations (14) and (16). 


R otor torque in shaft axes- The rotor aerodynamic torque equation (rof.l), 
which accounts for both acceleration and deceleration, is 
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wh ire 

S' 2 2 
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+ (0.00109 - 0.00270,^,5 - 3.13X* - 6.35X6, 75 - 1.930* 75 )n* 

- 0.133X0, + (-0.976X* - 6.38X0, ,, - 5.260*^,5)y'* (33) 


The aerocynamic t " que acting on the main rotor Qg„,, is calculated using main 
rotor parameters in (32) and (33). The torque applied to the fuselage by the 
riin rotor i function of Qnm and is determined by the engine and governor 
model: 


N, u " Q 

hub, 8 8 


(34) 


r . the main rotor, Og is equal to the engine torque Qeng* The tail rotor 
torque Qat* minted using tall rotor parameters in (32) and (33), is 
assumed to act directly on the fuselage so that Qg is equal to 


10 


Rotor hub aoaente to shaft axos- The hub momunts du« to flapping angle 
offsets are calculatad in shaft axea according to formulas obtained from 
reference 1. These formulas result from neglecting higher order terms In 
equations presented in reference 3: 


where 



(35) 


(36) 


are the flapping angles in shaft axes. 

Transformation of rotor moments to body axes- The rotor moments in shaft 
axes* given by equations (34) and (35), are transformed to body axes: 
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(37) 


The total moments applied to the fuselage by the rotor include the hub moments 
(37) and additional moments due to the location of the hub relative to the 
helicopter c.g.: 

■l" 

M 
N 

where the rotor forces are defined by equation (31). 



Tail rotor is hinge effect- The above model represents a rotor without 
a delta-three ( 63 ) hinge, such as the main rotor of a CH-53. However, the 
tall rotor has a 63 hinge, so that blade coning and flapping affect blade 
pitch; therefore, the model is modified accordingly. Asswning the changes in 
blade pitch due to flapping are small compared with those due to coning. 
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ot 


9^. “ « ► tan S.f. 
ct ot 


(39) 


where Oqj. is the value of collective pitch commanded by the control system. 
Note that the coning angle a^,, equation (21), is a function of 9 q; as a 
result, equations (21) and (39), for the tail rotor, should be solved 
simultaneously. 
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The Inputs and outputs of the rotor models are shown In figures 3^c) 
and 3(d). 


Engine and Governor Hodel 


An engine and governor model is included to provide a realistic time 
delay between aerodynamic rotor torque and the resulting reaction torque 
applied to the fuselage. This model was adapted from one used by Boeing Vertol 
(ref. 1); although it is not a model of a CH-S3 engine, it does provide the 
desired effects. 


This model, which includes the effects of a gas turbine, a power turbine, 
rotor Inertia, and shaft compliance, uses the reference rotor speed and 

the main rotor aerodynamic torque (eq. ( 32 )) to calculate the angular 

velocities of the main and tail rotors and the engine torque. Note that the 
engine torque Qeng equal to the main rotor shaft torque Qs in equa> 
tion (34). 
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(40) 


The constants and K<j represent the main rotor shaft compliance and damp- 
ing, respectively; note that the latter is required for computational stabil- 
ity. The term in the Qgen differential equation allows the model to 

hold reason.. ,)iy constant rotor speed under widely varying aerodynamic torques 
(ref. 1). 


The inputs and outputs of the engine model are shown in figure 3(e). 


Control System 

The control system model, which Includes the effects of pilot Inputs, 
control cross coupling, an automatic flight concrol system (AFCS) , and 
servo actuators, defines the main rotor collective pitch 0om» longitudinal 
and lateral cyclic pitch Bi and Aj, and tail rotor collective pitch command 
Occ* This model was obtained from an unpublished Sikorsky Aircraft report. 
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Pilot controls- 


where 
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(41) 


X' « 0 if X , < 2.54 cm 
col col 

or X . < 1.0 In. 
col 

and the term in parentheses Is limited to the range of -0.0349 rad (-2.0“) to 
40.419 rad (424. 0“). 


The pilot Inputs, in equation (41), are the displacements of the pilot 
controls relative to a nominal position. These positions are shown in the 
control rigging diagrams, figures 13(a)-13(d), as the zero displacement 
positions. The force characteristics of the pilot controls are given In 
table 2. 


AFCS inputs- The following features of the AFCS are not implemerted, 
directly, due to hardware limitations of the fixed-base simulator used in con- 
junction with this model: 


1. Trim adjustments for various c.g. locations 

2. Indicator of AFCS authority used 

3. Supplemental controller which changes effective collective stick 
position 


4. "Open- loop" pedal spring 

5. Lateral cyclic "stick pusher" 
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Effects of these features that are considered critical to the anticipated fly> 
ing tasks are included by modifying the AFCS model obtained from Sikorsky 
Aircraft. 

The absence of trim adjustments is compensated for by placing the c.g. at 
fuselage station 332 so that the mathematical model trims straight and level 
at 90 knots » with minimal AJPCS contribution to the longitudinal cyclic pitch; 
that is Bi^oklems caused by the lack of information on AFCS 

authority used are alleviated by (1) increasing the limits on the AFCS contri- 
bution to the tail rotor pitch command ^tafcs C^^Ble S); (2) removing &tafcs 
from the bracketed term in equation (41)* which is limited; and (3) selection 
of the c.g. fuselage station discussed above. The effects of the collective 
‘'.tick supplonental controller are not considered critical and. therefore, are 
not included. The "open- loop" pedal spring is represented by the integral 
term Kij/e the 6tafcs equation (42). The basic Effects of the lateral 
cyclic "stick pusher" are to provide the pilot with a stick force proportional 
to the deviation of roll attitude from its trim value, and to return the 
vehicle to its trim roll attitude when the pilot releases the stick. Since a 
"control loader" is not available in the fixed-base simulator* implementation 
of these effects required several changes to the AFCS model; the changes are 
described in detail below. 

In the original AFCS the roll trim ref<>rence is removed from the lateral 
channel, in equation (42), when the pilot places his feet on the 

pedals (activating a pressure sensitive .wlcch) prior to a lateral maneuver. 

If the pilot releases lateral stick pressure during the maneuver and keeps his 
feet on the pedals, the "stick pusher" moves the stick so as to regain the 
roll trim reference attitude. This characteristic is obtained by removing the 
roll trim reference from in equation (42) only when the lateral stick 

is displaced 1.27 cm (O.S in.) or more from its zero force trim position. 

Thus, the roll reference is removed when the pilot, by displacing the stick 
laterally, indicates a desire to maneuver; the reference is regained when the 
pilot releases his control force, allowing the stick to return to its zero 
force position. The control forces provided by the "stick pusher" during the 
maneuver are obtained by adding a bias proportional to the roll deviation from 
trim to Che lateral stick displacement (see eq. (42)). This causes the 

steady-state roll attitude deviation from trim to be proportional to lateral 
stick displacement from the zero force trim position and, therefore, propor- 
tional to the control force required by the pilot. The bias gain, Kjt,, in 
equation (42) corresponds to 0.14 N (0.08 lb) of pilot force per degree change 
in roll attitude. 

The control inputs from the modified AFCS model which contains altitude 
hold, heading hold, and turn coordination modes, are 
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i 



- <» 


altitude hold 


»«tc. ■ ] |- * *.»] 

. r fade In/out 1 F Kis „ . / fade in/out \„ ^ a\ 1 

xafca jclrcult no. 2 J^Tj8 + 1 


heading hold 


+ - h) 

altitude hold 


tafcs 


-(> 




fade In/out 
circuit no. 4 



T^S +1 h 


heading hold 


XI 


Hat 


lat 


[- 


/ fade in/out \ 
\circult no. 3/ 


turn coor- .ituude 
j dination_ 


. - h) 
hold 


*^*'*^trim^^trim “ *^^h 


(42) 


The fade In/out circuits are Intended to minimize the Introduction of 
large transients to the flight control system due to changes In the AFCS oper- 
ating mode. The gain of these circuits varies between zero and unity, accord- 
ing to the transfer functions listed In table 3. It should be noted that 
these transfer functions are only used to determine gain values, and do not 
represent actual filters. 

AFCS modes- The operational modes of the AFCS are controlled by the fol- 
lowing variables, which appear In equation (42). 


1^. > 1 Altitude hold mode engaged 

^ *0 Altitude hold mode disengaged 

1 . ■ 1 AFCS engaged 

* - 0 AFCS disengaged 

1 j 1 Pilot’s feet off pedals 

• 0 Pilot's feet on pedals 


trim 


1 

0 


Cyclic trim button released 
Cyclic trim button depressed 


I i 


} . 
I * 

i 

- s 


I . 
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1 


tc 


1 Above 60 knots and pilot's feet on pedals 

0 At or below 60 knots or pilot's feet off pedals 


^Xiat 


- 1 
• 0 


Lateral stick within 0.5 in. of zero force trim position 
Lateral stick beyond 0.5 in. of zero force trim position 


The values of the roll and heading trim angles, <^trlm "^trim respectively, 
are determined as follows: i^trim equal to the current ^he 

cyclic trim button is released; tpcrim is set equal to the current when 

the pilots feet move off the pedals. 

The AFCS and the altitude-hold mode are activated by switches on the 
Instrument panel. The heading-hold and turn-coordination modes are controlled 
by airspeed and location of the pilot's feet (either on or off the pedals). 

The I ading-holc mode is engaged whenever the pilot's feet are off the pedals, 
regardless of airspeed. The turn-coordination mode is engaged only when the 
pilot's feet are on the pedals and the airspeed is greater than 60 knots. The 
operation of these modes is summarized in table 4. 


AFCS authority limits- The authority of the AFCS is limited so that it 
can be overridden by the pilot. This is accomplished by limiting the control 
inputs from the AFCS, equation (42), to the values shown in table 5. In the 
expressions for and limits are imposed prior to the addi- 

tion of the altitude-hold terms. 


Servo actuators- The primary servo actuators transform the main rotor 
control commands, given in equation (41), into swashplate angles and blade 
collective pitch. The following model of these servos was obtained from 
Sikorsky Aircraft. 


"9' . 
om| 

b; 


2 “toS 

0 ) e ° 
n 


e 


A! 


(s* + 2?W^S + 'ilj^^)(TS + 1) 


-•( 8 ) 


om 
A, 


(43) 


L."l 


( 8 ) 


A model of the tail rotor servo was not obtained from Sikorsky; therefore, it 
was assumed that 


Approximations for real-time simulation- During use of this helicopter 
model in real-time guidance and navigation studies it may be desirable to 
neglect some of the relatively high-frequency dynamics — specifically, the 
relatively small time constants, and Xj, in equation (42), and the servo 
dynamics, equation (43). 

The inputs and outputs of Che control system model are shown in fig- 
ure 3(f). 
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Equations of Motion 


The helicopter equations of motion are given in body axes with respect to 
a flat, nonrotating Earth. The helicopter is considered a rigid body with 
mass symmetry about the xj, - S}) plane. The effects due to the engine angular 
momentum are neglected. 


Translational acceleration- The translational equations of motion are 


"'h/e 
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“x" 
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r, h 

Lwj 

eg, h 

P 0- 
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cg» h^ 


where 


n/e 


cos 6 cos l|» 

sin sin 3 cos <|> 
• CCS (|i sin 

cos $ sin 6 cos i|i 
L + sin sin 


cos 0 sin 

cos cos ’ll 

■f sin 41 sin 6 sin 

cos 4 > sin 6 sin tp 

- sin 4 cos 


-sin 9 
sin 41 cos 6 

cos 4 > cos 6 


(A5) 


(46) 


and 4h> ^h* 'I'h Euler angles that define Che orientation of the 

body axis system (fig. 3). The fuselage aerodynamic forces are given by 
equation (10), and the rotor forces, which Include those due to both main and 
tail rotors, are given by equation (31). Equation (45) can be rearranged to 
yield 



m 






(47) 


Inertial velocity and position- The inertial velocity, in body coordi- 
nates, is obtained by integrating equation (47), with respect to time, subject 
to appropriate initial conditions. The inertial velocity in Earth axes is 


V 


V 




V 


V 

V 



w m 

eg, e 



'eg, h 


(48) 


The position of the helicopter, in Earth coordinates, is determined by inte- 
grating equation (48) with the appropriate initial conditions. 
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( 49 ) 



Angular acoleration- The rotational equations of motion are 


where 



(50) 


The fuselage aerod aamlc momenta are given by equation (11)* and the rotor 
moments • which Include chose due to both the main and tall rotors » are given 
by equation (38). Equation (50) can be rearranged to yield 
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f. h 
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-r 
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0 

-P 
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P 
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r^ 
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(51) 


Angular velocity and orientation- The angular velocity, in body axes, is 
obtained by integrating equation (51), with respect to time, subject to thr 
appropriate Initial conditions. 


The helicopter Euler angles are determined by integrating 


where 


’V 


mm m 

P 

0 

■ 

o 

1 

>-* 

q 



JC, 



1 

0 

-sin 6 

“ 

c - 

0 

cos 

sin ^ cos 
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i) 

-sin 

cos cos 
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(52) 


( 53 ) 


The inputs and outputs of the equations of motion are shown in fig- 
ure 3(g). 
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MODEL VALIDATION 


Th« mathtmatlcal model it validated by comparing its response to chat of 
an actual helicopter and by a qualitative comparison of the handling charac- 
teristics made by experienced pilots. 


Time History Comparisons 

The most readily available flight data were from an HH-53C, an Air Force 
version of a CH-SSC* which has two external fuel tanks. Since the HH-S3C 
response time histories given In reference 6 were obtained with these tanks 
fully the helicopter Inertias In the model were modified accordingly In order 
to provide a more realistic comparison of responses. The following modified 
parameters were calculated using data supplied by Sik. rsky Aircraft. 

Ixxyh • 56,367 kg-m^ 

■ 41,553 slug-ft^ 

lyy,h - 268,709 kg-m^ 

- 198,090 slug-fi.^ 

l2Z,h « 248,745 kg-m* 

- 183,373 slug-ft^ 

1X2, h ■ 28,400 kg-ffl^ 

- 20,936 slug-ft^ 

The aerodynamics effects of the external tanks were not known and, there- 
fore, not incorporated into the model, leference 6 contained HH-53C time 
histories for pulse type inputs to the longitudinal cyclic, lateral cyclic, 
and pedals, at the following flight condition with the AFCS both on and off. 

Air speed ■ 113 knots 

Altitude • 7C00 ft 
Main rotor speed * 185 rpm 
Gross weight - 41 ,000 lb 
FSCG - 328 

Atmospheric temperature - -18’ C 

The response time history cf the CH-53 model at the above flight condi- 
tion was obtained using a "dynamic check" routine. This routine provided the 
model with flight control inputs that approximated those of the HH-53C. Also, 
this routine was used to control the operational modes of the AFCS, as will be 
discussed later. The time histories of the CH-53 model and the HH-53C are 
compared with the AFCS on; this is done because the model will normally be 
operated in this mode for terminal-area studies. 

A comparison of the .espouses to a forward longitudinal cyclic pulse, 
shown in figure 14(a), Indicates good agreement for the Euler angles and fur 
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the body-axis angular velocities. This is also th'C casa for the responses to 
a right lateral eye, lie pulse* shown in figure 14(b). Here -it is assumed that 
the HH-S3C response was. obtained with the cyclic trim button depressed* since 
the roll attitude does not return to aero after the pulse.. This condition was 
simulated in the model by using the dynamic check routine to set Itrim * ^ 
in equations (42). The responses to a right. ped^^l pulse did not compere as 
well as those for the previous inputs.*- The model produced much larger atti- 
tude excursions than indicated for tha l}H-b3C -response. A reasonable compari- 
son* shown in figure 14(c), was obtained by raising the damping gain K 19 
-from 0.57'3 to 1.50. 


r' 




Possible sources of, the discrepancy are the unmqdeled aerodynamics of the 
HH-53C external fuel tanks* and features of the AFCS which were modified or 
not Included due to limitations of the fixed-baSe simulator. An attempt jrfas 
made to compensate for differences between the actual and modeled AFCS by coj»- 
trolling the mode"' of the latter with the dynamic check routine. For the 
response to„ a peual pulse*. this routine simulated the AFCS transformation from 
the heading-hold mode, -to the tumrcoordination mode by setting Ip^^j >= 0 and 
Itrim ® ® equations (42), Later review indicated hat this is a poor 
method for simulating the mode transformation. The above method completely 
removes the roll trim reference when the pedal pulse is initiated. Actually, 
this re'erence should fade out with a l-sec time constant and* therefore, a 
rabre realistic simulation would keep Itrim “ ^ ^xlat' ® ^ 

Iped ® 0 at the beginning of the pulse. It may also be desirable to eliminate 
effects of the "stick pusher" by setting K 24 =0 in equation (42). Because 
the method used in simulating the AFCS mode transformation served to prema- 
turely remove an attitude error signal* it probably increased the attitude 
excursions of the model and* therefore* may have contributed to the response 
discrepancy. 


Pilot Comments 

A qualitative evaluation of the mathematical model was made by two pilots 
using a fixed-base simulator with visual scene. These evaluations were to be<- 
made ccmsidering the intended use of the model* that is* terminal-area guidance 
and navigation studies. 

The control forces and general feel of the flight controls were satisfac- 
tory, although' the absence of breakout and gradient forces* with the cyclic 
trim button depressed* results in stick-jump and a tendency to overcontrol. 

The absence of cyclic beeper trim and collective and pedal parallel servos did 
not degrade the model for its Intend.-d use. The basic AFCS functions were 
primarily evaluated in forward flight at approach speeds (90-120 KIAS) . The 
retentit i of trimmed airspeed and pitch and roll attitude was excellent. The 
AFCS modes* altitude-hold* heading-hc'’d* and furn-coordiiiation operated satis- 
factorily. Also* the AFCS modifications* made :o include effects of a lateral 
"stick pusher," provided responses that were much more ch.nracteristic of the 
aircraft. Although* not required* low-speed flight and hover were also eval- 
uated. .The attitude- and heading-hold features operated very well during 
decelerating approaches to a 50-ft hover. The collective increases and nosv'- 
high attitudes required during deceleration were similar to those of the CH-33 



aircraft. Abova 10 knota, turna ware aaally coordinated with tha padala. At 
lower apoadBi in forward and aidaward flight and in hover, praciao haading and 
hover control raqulrad much cloaer pilot attantion to tha turn coordinator. 
Thla waft mainly dua to inaufficiant motion cuaa from tha viaual acana. 

It waa concluded that the flying qualltiaa of thla modal wara qualita'- 
tivaly repreaentativa of the actual aircraft, within the limitation of a 
fixed -baaa aimulator. 


CONCUSSIONS 


The mathematical model of a CH'-53 helicopter dvacrlbed In thla report waa 
developed for real-tlmo piloted aimulatlon in support of terminal-^area guidance 
and navigation studies. This model is based on modified nonlinear classical 
rotor theory and nonlinear fuselage aerodynamics, limitations of the fixed* 
base simulation facility prevented direct implemuntation of several features 
of the automatic flight control system (AFCS). Tha eff ^cts of these features 
considered critical to the anticipated flying tasks were included by modifying 
the AFCS equations. The model was validated by comparing its response with 
actual flight data and by a qualitative comparison of the handling character- 
istics made by experienced pilots. These comparisons indicated that the model 
Is sntisfuctory for terminal-area guidance and navigation stvidles. 
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TABLE 2.~ FORCE CHARACTERISTICS OF PILOT CONTROLS 


Pilot control 

Force 1 

Breakout » N (lb) 

Gradient. N/cm (Ib/ine) 

Longitudinal cyclic 

■nnRnRn 

2.3 (1.3) 

Lateral cyclic 


1.4 (0.8) 

Collective 


0 (0) 

Pedals 


5.3 (3.0) 


TABLE 3.- GAINS OF FADE IN /OUT 


1 Circuit number I 

1 

2 

3 

4 

^afcs 

^afcs 

^xiat 

^afcs 

tjS + 1 

TgS + 1 

TyS + 1 

T,S + 1 


TABLE 4.- SUMMARY OF HEADING-HOLD AND TUSN-COORDINATION MODE LOGIC 


1 Operating condition 

Feedback 

information 

Ifode 

\WMS3S^3KtM 

Pilots feet 


Heading 

Roll rate 

Lat. acc. 

Above 60 knots 

On pedals 

Off 

Off 

On 

On 

Tum-coord . 


Off pedals 

On 

On 

Off 

Off 

Head* -hold 

Below 60 knots 

On pedals 

Off 

Off 

Off 

Off 



Off pedals 

On 

On 

Off 

Off 

Head. -hold 


TABLE 5.- AUTHORITY LIMITS OF AFCS 


AFCS input 

Limits, rad ^deg) 

^mafcs 

±0.0227 rad, (±1.3) 

®iafcs 

±0.0454, (±2.6) 

^lafcs 

±0.0209, (±1.2) 

^tafcs 

±0.0332, (±1.9)^^ 


^hls limit Increased to ±7° 
in the simulation. 





























Sikoriikv Cn-5i 


Figure 





Figure 2.- Helicopter body axes, shaft axes, control axes, and flapping 

angle definitions. 
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(a) Helicopter simulation. 


Block diagram of helicopter model and input-output diagrams of 
individual component models. 
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(d) Tail Rotor (subscript t) model. 


Figure 3.- Continued. 
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(e) Engine and "overnor model. 


Figure 3.- Continued. 
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(b) Longitudinal cyclic 


Figure 13.- Continued. 








LATERAL CYCLIC PITCH 
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Figure 14.- Flight test-math model comparisons. 
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